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An extension of the single-freeze-out model with thermal and geometric parameters dependent 
on the spatial rapidity, atu, is used to describe the rapidity and transverse-momentum spectra of 
pions, kaons, protons, and antiprotons measured at RHIC at v / sjvjv = 200 GeV by the BRAHMS 
collaboration. THERMINAT0R is used to perform the necessary simulation, which includes all resonance 
decays. The result of the fit to the rapidity spectra in the range of the BRAHMS data is the expected 
growth of the baryon and strange chemical potentials with the magnitude of atu, while the freeze- 
out temperature is kept fixed. The value of the baryon chemical potential at qm ~ 3, which is the 
relevant region for particles detected at the BRAHMS forward rapidity y ~ 3, is about 200 GeV, 
i. e. lies in the range of the values obtained for the highest SPS energy. The chosen geometry of the 
fireball has a decreasing transverse size as the magnitude of an is increased, which also corresponds 
to decreasing transverse flow. This feature is verified by reproducing the transverse momentum 
spectra of pions and kaons at various rapidities. The strange chemical potential obtained from the 
fit to the K + J K~ ratio is such that the local strangeness density in the fireball is compatible with 
zero. The resulting rapidity spectra of net protons are described qualitatively in the model. As 
a result of the study, the knowledge of the "topography" of the fireball is achieved, making other 
calculations possible. As an example, we give predictions for the rapidity spectra of hyperons. 
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INTRODUCTION 



The study of particle abundances has been a major 
source of information concerning heavy- ion collisions. In 
fact, the agreement of the particle ratios with simple 
predictions of statistical models is a key argument for 
early thermalization of the formed system [l|, 0, 0, 0]. 
Up to now the numerous studies of the particle ratios 
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[2U [22l . [23| were falling into two basic categories: the 
so-called Air studies at low energies (SIS, AGS) and the 
studies at mid-rapidity for approximately boost-invariant 
systems at highest energies (RHIC). The An studies in- 
volve three-momentum integrals of the statistical distri- 
bution functions, with the multiplicity of species i given 
by Ni = V J d?pfi{y/ mf + p 2 ), thus providing informa- 
tion on volume-averaged thermal parameters of the sys- 
tem in a very simple way. The inclusion of resonance de- 
cays 0, HE HI] , crucial for the success of the approach, 
is also straightforward, since the detection of the prod- 
ucts with full angular coverage is insensitive to the decay 
kinematics or flow effects; once the resonance has de- 
cayed, its products are registered. The other simple sit- 
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uation arises when the system is nearly boost-invariant. 
To a sufficiently good accuracy this is the case at mid- 
rapidity for the RHIC energies, where the particle yields 
dN/dy change only by a few percent in the rapidity win- 
dow \y\ < 1. The assumption of the boost-invariance of 
the fireball leads again to very simple formulas. Although 
the particles detected at mid-rapidity are collected from 
various parts of the fireball, not only from the very cen- 
tral region, their ratios are the same as in the 47r calcula- 
tion. This is because for dNi/dy = J d 2 p±d 3 Ni/(d 2 p±dy) 
we have from the boost invariance 12711 



dNi/dy = J dy dNi/dy = Nj_ 
dNj/dy ~ JdydNj/dy ~ Nj 



(1) 



This obvious general formula finds an explicit manifesta- 
tion in specific boost-invariant models. The result also 
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holds when resonance decays are included, see Ref. 28 1 
for a derivation in the framework of the Cooper- Frye 29 1 
formalism. 

When the system is not boost invariant the above sim- 
plifications no longer hold. The situation is illustrated 
in Fig. [T] Particles detected at a given pseudorapidity 77 
(parallel lines in the figure) originate from different pieces 
of the fireball (gray blobs). Thermal conditions (tem- 
perature, chemical potentials, flow) change from piece 
to piece, which must be properly included. In addition, 
the effects of the longitudinal flow (indicated by arrows) 
must be incorporated, and the kinematics of resonance 
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FIG. 1: (Color online) Emission of particles from a boost- 
non-invariant fireball. The horizontal (vertical) axis indicates 
longitudinal coordinate z (transverse coordinate p) . Particles 
emitted with the same value of pseudorapidity r\ originate 
from different regions of the fireball indicated by the gray 
blobs (they also are emitted at different times). The thermal 
conditions and flow (indicated by arrows) in these regions 
are different. The solid lines indicate tracks of primordial 
particles, while the dashed lines show products of resonance 
decays. The dashed ellipse indicates the relevant region for 
a give 77, which spans about two units of the spatial rapidity 



which is the relevant region for particles detected at the 
BRAHMS forward rapidity, y ~ 3, is about 200 GeV. 
This value is in the range of the values of the thermal fits 
for the highest SPS energy, thus we confirm the recent 
findings by Roehrich 31( that the thermal conditions at 
RHIC at forward rapidities, y ~ 3, correspond to the 
SPS conditions at mid-rapidity. Details of our procedure 
of determining the dependence of thermal parameters on 
a || are explained in Sect. Our strategy, as usual, is 
to fix the features of the fireball with the well-measured 
spectra of particles: pions, kaons, protons and antipro- 
tons. The strange chemical potential obtained from the 
fit to the K + /K~ ratio is such that the local strangeness 
density on the freeze-out hypersurface E is compatible 
with zero. The experimental [32| rapidity spectra of net 
protons, p — p, are reproduced qualitatively in the model, 
displaying the correct shape but overshooting the data 
at larger rapidities by about 50%. The chosen geome- 
try of the fireball incorporates a decreasing transverse 
size as |aii| is increased, which simultaneously results in 
a decreasing transverse flow. This choice is verified in 
Sect. IIIII by reproducing the spectra d 2 N/(2irpT<lpTdy) 
at a fixed y of pions and kaons at ^/sJTn = 200 GeV 
from the BRAHMS collaboration [33] . These transverse- 
momentum spectra exhibit slopes which become steeper 
with rapidity. 

As a result of our study, we obtain the "topography" 
of the fireball, which can be the ground for other more 
detailed studies, discussed in the Conclusion. 



decays (dashed lines) becomes relevant. Although the 
resulting formalism for particle spectra remains concep- 
tually simple and is based on the standard Cooper-Frye 
treatment, the calculation is no longer semi-analytic and 
a full-fledged simulation is necessary to accomplish the 
goal. 

In the analysis of this paper we use THERMINATOR - the 
THERMal heavy IoN generATOR [H, to generate the 
Monte Carlo events in a suitably modified single-freeze- 
out model of Ref. [27| . The extension to the boost-non- 
invariant case consists of two basic elements. The first 
one (geometric) is the choice of the shape of the freeze-out 
hypersurface E and collective expansion. The other one 
incorporates the dependence of the thermal parameters 
on the position within the hypersurface E. Specifically, in 
our treatment the transverse size and the chemical poten- 
tials depend on the spatial rapidity an = arctanh(z/t), 
where z and t are the longitudinal and time coordinates 
on the freeze-out hypersurface. Although the boost-non- 
invariant model has quite a few parameters, as listed 
at the end of Sect. [IIJ they can be fitted independently 
to various combinations of the data, leaving little free- 
dom. For instance, the an dependence of the baryon 
and strange chemical potentials is fixed with the ratios 
of protons to antiprotons and K + to K~ . The result in 
the range of the BRAHMS data is the expected growth 
of the baryon and strange chemical potentials with | 1 1 1 - 
The value of the baryon chemical potential [Ib at an ~ 3, 



II. THE SINGLE FREEZE-OUT MODEL 

The single- freeze-out model is described in detail in 
Refs. [13, HI, . Here we review the main assumptions 
and the formalism of describing the expansion and par- 
ticle decays. 

1. At a certain stage of evolution of the fireball the 
thermal equilibrium between hadrons is reached. 
Most probably, hadrons are "born" already in such 
an equilibrated state. The local particle phase- 
space densities have the form of the Fermi-Dirac or 
Bose-Einstein statistical distributions. The parti- 
cles generated at freeze-out are termed primordial. 
For simplicity of the model we do not include the 7 
non-equilibrium factors of Ref. [35| used in several 
recent analyses [3, USUI. 

2. The thermodynamic parameters are the freeze- 
out temperature T and three chemical potentials: 
baryon, /is, strange, /15, and /i/ 3 , related to 
the third component of isospin. In a boost-non- 
invariant model the values of these parameters de- 
pend on the position within the freeze-out hyper- 
surface E. 

3. In the boost-non-invariant model the shape of the 
fireball is nontrivial in the longitudinal direction. 
In this paper we retain the azimuthal symmetry. 
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4. The velocity field of the collective ex pan sion is cho- 
sen in the form of the Hubble flow [3a ] , providing 
the longitudinal and transverse flow to the system. 
Again, in the boost-non-invariant model the func- 
tional form of the velocity field may depend on the 
longitudinal position. 

5. The evolution after freeze-out includes decays of 
resonances which may proceed in cascades. All res- 
onances from the Particle Data Tables [37| are in- 
corporated. 

6. Elastic rescattering among particles after the chem- 
ical freeze-out is ignored and the model may be 
viewed as an approximation to a more detailed evo- 
lution, taking into account different time scales for 
various hadronic processes (see [38| and references 
therein) . 

The single freeze-out concept complies to the explo- 
sive scenario at RHIC fl4j . Moreover, the approach 
reproduces very efficiently the particle abundances, the 
transverse-momentum spectra, including particles with 
strangeness [HI , produces very reasonable results for the 
resonance production [3i|, the charge balance functions 
in rapidity 0, the elliptic flow SI, the HBT radii 0, 
and the transverse energy [HI, S3 SH- Recently with 
the help of RQMD it was found (cf. Fig. f 6 and 17 
of Ref. [47J) that the elastic rescattering effects are not 
significant for the mid-rapidity pt spectra of pions. One 
can understand this as follows: resonance decays "cool" 
the spectra [13] • As the result, the original ~ 165 MeV 
spectra from the chemical freeze-out look, after feeding 
from resonances, approximately as ~ 130 MeV spectra, 
which would be obtained at the the lower temperature 
of the thermal freeze-out. Thus clastic rescattering be- 
comes innocuous. Certainly, more studies are desirable 
here, in particular an "afterburner" performing elastic 
rescattering could be run on top of our simulation. It 
would help to achieve a more accurate collision picture, 
with the elastic rescattering processes taken fully into 
account. 

Popular choices of the freeze-out hypersurface and the 
collective velocity field are discussed in detail in Ref. [48[ . 
In this work we modify in a very simple way the orig inal 
boost-invariant single- freeze-out model of Ref. [27l.[28l.l33] 
by making the transverse size of the fireball dependent on 
the spatial rapidity. We use the freeze-out hypersurface 
parameterized as 
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(2) 



The parameter a\\ is the spatial rapidity, while a±_ is re- 
lated to the transverse radius as 



The four- velocity field is chosen to follow the Hubble law 

u ' 1 =x^/t. (4) 

We note that the longitudinal flow is v z = tanhan = z/t 
(as in the one-dimensional Bjorken model 49]), while the 
transverse flow (at z = 0) has the form v p — tanhaj_. 

The new element of the parameterization of this paper, 
which implements the departure from boost invariance, 
is the selection of boundaries for the fireball. In the boost 
invariant model p was limited by the space-independent 
parameter p max , or < a± < a™ ax . Now we take 



< ol < a]™ x 
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The interpretation of this formula is clear: as we depart 
from the center by increasing |a|||, we simultaneously re- 
duce a± 1 or p max . The rate of this reduction is controlled 
by a new model parameter, A. Since in our model the 
flow is linked to the position via Eq. ([JJ, we also have 
less transverse flow as we increase |aii|. This feature will 
show in the px spectra presented in Sect. IIIII We may 
also use more conveniently the parameter 
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(0). 



(6) 



Thus the geometry and expansion of the fireball is de- 
scribed by three parameters: t, pmix, and A. 

With the standard parameterization of the particle 
four-momentum in terms of rapidity y and the transverse 
mass mi, 

p^ = (m±coshy, p± cos ip, p± sin tp, m±smhy) , (7) 
we find with Eqs. © and @ 

p-u = m^cosh(a_L)cosh(a|| ~y) — _p_Lsinh(aj_) cos(</>— ip), 

(8) 

and 

d 3 T, ■ p — dctud(f> pdp x 



i±V t2 + P 2 cosh(a|| - y) - p±p cos(0 - (p) 
T i da^d(f)sinha±dct±p ■ u (9) 



where d 3 E^ is the volume element of the hypersurface. 
With the assumed azimuthal symmetry the Cooper-Frye 
[291 ] formalism then yields the following expression for 
the momentum density of a given species of primordial 
particles: 



d 2 N 
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Tsinha 
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2npxdpT dy 



r Q T ax ( Q n) r 27r 
dan / da± / d<f> 
> Jo Jo 

xp ■ uf (ftp ■ u - (3{j,(a\\)) , (10) 

1 1 



(27r) 3 exp z ± 1 ' 
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where p ■ u from Eq. @ is taken at ip = 0, /(z) is the 
statistical distribution function (with + for fcrmions and 
— for bosons), = 1/T, and 

(i(a\\) = Bfi B (at\\) + Sfis(oc\\) + I 3 fii a (a\\), (11) 

with B, S, and /3 denoting the baryon number, 
strangeness, and the third component of isospin of the 
particle. Thus we admit the dependence of chemical po- 
tentials on the spatial rapidity. This is of course neces- 
sary if we wish to describe in the framework of a sta- 
tistical model the increasing density of net protons as 
we move from mid-rapidity towards the fragmentation 
region. 

The temperature T also in general depends on an. 
The best model-building strategy here would be to 
use the universal Cleymans-Redlich chemical freeze-out 
curve fl2| in the \±b-T space (for a recent status see 
Ref. jl8l. 50|). That way the functional dependence of 
\xb on a || induces unambiguously the dependence of T 
on an. In this work, however, we apply the model for 
not too large values of the rapidity, \y\ < 3.3, and it will 
turn out that the obtained values for \xb are less than 
~ 250 MeV. The universal freeze-out curve gives from 
fi b = to fi B = 250 MeV a practically constant value of 
T. For instance, at SPS {Pb + Pb, TsTTn = 17 GeV) one 
has [IB] T = 164 MeV, /jb — 229 MeV, us = 54 MeV, 
and fii 3 = —7 MeV, with the value of T equal within 
errors to the RHIC value of 165 MeV. For this reason in 
the present analysis we fix the freeze-out temperature at 
a constant (independent of the spatial rapidity) value, 



T = 165 MeV. 



(12) 



If modeling were made for larger values of the rapidity 
and/or lower collision energies, the dependence of T on 
qt|| should be incorporated according to the prescription 
based on the universal freeze-out curve. Eventually, we 
expect that when the fragmentation region is approached, 
T becomes very small and \xb reaches the value of the 
order of 1 GeV. 

Another qualitative argument for the approximate con- 
stancy of T at moderate values of |a||| may be inferred 
directly from the BRAHMS data. From the measured ra- 
pidity spectra (cf. Fig. [5]), obviously, the yields of pions 
and kaons decrease with y. Thus, one needs to decrease 
the size of the emitting source, decrease the freeze-out 
temperature, or both. The decrease of temperature af- 
fects more strongly the particles with larger masses, since 
the thermal factor is approximately exp(— m 2 +p 2 /T). 
Therefore, if we introduced variation of the temperature 
with the spatial rapidity, it would result in a faster drop 
with y of the pion yields compared to the kaons. Cer- 
tainly the data excludes this situation, since the ratio of 
dNir/dy to dNx/dy is within a few percent independent 
of y in the BRAHMS rapidity coverage. Therefore we 
must keep T constant (within a few percent), and the re- 
maining possibility is the decrease of the source size with 
| q; |||- Resonance decays complicate the above qualitative 
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FIG. 2: (Color online) Top: the model baryon and strange 
chemical potentials plotted as functions of the spatial rapid- 
ity. Parameters of Eq . (|13p are obtained from the fit to the 
BRAHMS data [32l |33|. The points represent a naive calcu- 
lation based of Eq. (I15|l . Bottom: the ratio of the baryon to 
strange chemical potentials, /xb/ms- 



argument, but with the help of a numerical simulation we 
confirm it. Another way of providing the drop of yields 
with rapidity is to incorporate the 7 non-equilibrium fac- 
tors [35j dependent on an, which may dilute the system 
as |ay | increases. We do not explore this possibility here. 

For convenience, we parameterize functionally the de- 
pendence of the chemical potentials at low values of \a\\ | 
as follows: 



Ati(aii) = Hi{0) 



1 + Act 



2.4 



B,S,I 3 



(13) 



The chosen power of 2.4 works somewhat better than 2. 
Of course, any convenient and sufficiently rich parametric 
form is admissible here, as it is fitted to the data (see 
Sect. IIII[) and the introduced parameters effectively are 
not free. By "low" |aii| we mean the values relevant to 
the BRAHMS data, covering \y\ < 3.3. 
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Formula provides the spectra of the primordial 
particles. The following evolution of the system con- 
sists of free streaming, with resonances decaying into the 
daughter particles. We use THERMINATDR to perform the 
simulation. The code incorporates all the **** and *** 
resonances. Following the scheme of SHARE 51] it ex- 
cludes all * resonances, and practically all ** resonances 
listed in the Particle Data Tables [37|. Each resonance 
decays at the time controlled by its lifetime, 1/r. In the 
resonance's rest frame the decay at time t occurs with the 
probability density rexp(— Ft). The two-body or three- 
body decay channels are incorporated and the values of 
the branching ratios are taken from the Particle Data Ta- 
bles. Heavy resonances may of course decay in cascades. 

Let us summarize the model parameters. We have the 
universal freezeout temperature T, three chemical poten- 
tials at mid-rapidity, /is(0), /xs(0), /i/ 3 (0), three param- 
eters Ab , As , and Aj 3 of Eq. (fT5|) describing the depen- 
dence of the chemical potentials on the spatial rapidity, 
and three geometry /expansion parameters: the proper 

time r, the transverse size at mid-rapidity, pmix, and the 
parameter A controlling the spatial rapidity dependence 
of the transverse size. Except for T taken to have the 
value (|12p obtained in earlier thermal analyses of parti- 
cle ratios [27| , the remaining parameters are fitted to the 
BRAHMS data [Hill for the double differential spectra 
d 2 N/(2Trp T dp T dy). 



III. RESULTS 

We first describe our fitting strategy, which with many 
parameters present must be done with care. We wish 
to have a good starting point for the parameters de- 
scribing the chemical potentials. Practice shows that 
to a very good approximation the statistical distribu- 
tions are very well approximated by the Boltzmann fac- 
tors. Then the integrand of Eq. (|TTJ|) contains the fac- 
tor exp[— /3m_LC0sh(a:j_)sinh(a|| — y) + f3/j,(a\\)}. Because 
of this the relevant integration range in aii is sharply 
peaked around aii ~ y (the half-width is about one unit 
of a||) and the chemical potentials entering the formula 
are taken approximately at ^i(ot\\) ~ t*i(y)- Thus the fac- 
tors exp[/3 fi(y)} can be taken in front of the integration. 
If it were not for the resonance decays which modify the 
result (and are included in the full simulation) we would 
have to a good approximation the following relations: 



V 
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El 
K~- 

7T~ 



exp(2/3/i B ), 

exp(2/3/is), (approximate formulas) 
exp(2/3/i /3 ). (14) 
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The symbols on the left-hand side denote the ratios of 
yields of the specified particles at a fixed y and integrated 
over pt- These are known form the data, thus we can 



FIG. 3: (Color online) Dependence of particle ratios on rapid- 
ity: top - p/p, middle - K + /K~ , bottom - n + /tt~ . The open 
triangles are the Brahms data [32 . |33| , while the filled circles 
show the result of the model simulation with THERMINAT0R. 
The model parameters are given in Eq. (|17[l . The data for p 
and p are not corrected for the feed-down from weak-decays 

a. 



invert 



^b(jj) = -^TXogijp/p), (approximate formula) (15) 
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and so on. With the help of this form we set the starting 
values of the parameters Pi(0) and A iy which are then it- 
erated. The iteration proceeds as follows: for a given set 
of parameters we run the full THERMINATDR simulation, 
which generates events. We first optimize the baryon- 
number parameters /is(0) and As with the help of the 
ratio of the p and p rapidity spectra, then the strangeness 
parameters fis(0) and As using the K + to K~~ ratio, 
then we go back again to the baryon parameters, etc., 
and loop until a fixed is reached. The isospin parameters 
fij 3 (0) and Aj 3 are consistent with zero and thus irrele- 
vant. The A parameter is fixed with the pion rapidity 
spectra dN 7r ± j dy. The optimum value is 

A = 3.33. (16) 

The result of our optimization for the chemical poten- 
tials is shown in Fig. [2j The optimum parameters are: 

fi B (0) = 19 MeV, ns(0) = 4.8 MeV, /i/ 3 (0) = -1 MeV, 
A B = 0.65, A s = 0.70, A h = 0. (17) 

We observe the expected behavior for the baryon chem- 
ical potential, which increases with |aii|. The value at 
the origin is 19 MeV, somewhat lower than the ear- 
lier mid-rapidity fits made in boost-invariant models in 
Refs. [TEl \2l\. yielding 26 MeV. The lower value in our 
case is well understood. The previous mid-rapidity fits 
include the data in the range \y\ < 1. This range col- 
lects the particles emitted from the fireball at |a||| < 2, 
hence the value of \xb in the previous mid-rapidity fits is 
an average of our [Ib{oi\\) over the range, approximately, 
|aii| < 2, with some weight proportional to the parti- 
cle abundance. This qualitatively explains the effect of 
a lower value of our /xs(0) than in the boost-invariant 
models. A similar effect occurs for ps- We do not incor- 
porate corrections for the feed-down from weak decays 
(em i.e. all decays are included), since this is the policy 
of Ref. [32j for the treatment of p and p. 

We note that at a\\ =3 the value of ps is 200 MeV, 
more than 10 times larger than at the origin. This value 
is comparable to the highest-energy SPS fit (a/sat n = 
17 GeV), where p B ^ 230 MeV. The behavior of the 
strange chemical potential is qualitatively similar. It also 
increases with \a» |, growing form 5 MeV at the origin to 
50MeVata|| =3. The ratio pb{o-\\) / Psi a \\) is ver y close 
to a constant, ~ 4 — 3.5, as can be seen in the bottom 
panel of Fig. [2] 

The points in the top panel of Fig. [2] show the result of 
the naive calculation of Eq. (|14ll5p . We note that these 
points are very close (in particular for the strangeness 
case) to the result of the full-fledged fit of our model. 
This is of practical significance, since the application of 
Eq. (|14I15[) involves no effort, while the model calculation 
incorporating resonance decays, flow, etc., is costly. 

There is another important point. In thermal models 
one may obtain the local value of the strange chemical 
potential, ps, aA, & given ps with the condition of the van- 
ishing strangeness density, ps = 0. The result is shown in 
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FIG. 4: (Color online) Comparison of the strange chemical 
potential obtained from the fit to the data (solid line) and 
from the condition of zero local strangeness density, ps — 
(dashed line). 



Fig. |U where we compare the strange chemical potential 
obtained from the fit to the data (solid line) and from 
the condition of zero local strangeness density at a given 
Pb{&\\)- The two curves turn out to be virtually the 
same. This shows that the net strangeness density in our 
fireball is, within uncertainties of parameters, compatible 
with 0. This is not obvious from the outset, as the con- 
dition of zero strangeness density is not assumed in our 
fitting procedure. Although this feature is natural in par- 
ticle production mechanisms, in principle only the total 
strangeness, integrated over the whole fireball, must be 
initially zero. Variation of the strangeness density with 
ck|| is admissible, but turns out not to occur. 

Figure [3] shows the quality of our fit for the parameters 
of chemical potentials, Eq. fTTj) . We show the measured 
ratios oip/p, K + /K~ , and ir + /tt~ as a function of rapid- 
ity y [13, HH and the results of fit made with help of the 
simulation with THERMINAT0R. We note a very reasonable 
agreement. The error bars on the model points are sta- 
tistical errors due to the finite size of the sample (we use 
2500 simulated events in this plot). The flat character of 
the 7r + /7r~ ratio indicates that the value of the isospin 
chemical potential is consistent with zero at all spatial 
rapidity values. 

In Fig. [5] we show the comparison of obtained rapid- 
ity spectra of ir + , K + , and K~ to the experimental 
data. The experimental yields for the pions are corrected 
for the feed-down from the weak decays as described in 
Ref. [33j . For that reason for the case of tt + we give the 
model predictions with the full feeding from the weak 
decays (solid line) and with no feeding from the weak 
decays at all (dashed line). We note a quite good quanti- 
tative agreement, with the data falling between the two 
extreme cases. We recall that the behavior on rapidity of 
dN/dy is controlled by the A parameter of Eq. ([5]) . The 
spectra of 7r~ are not shown, since they are practically 
equal to the case of tt + . The spectra of K + and K~ are 
also quite well reproduced. 

Figure [5] displays the rapidity spectra of protons and 
antiprotons, as well as their difference p — p, i.e. the 
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FIG. 5: (Color online) Rapidity spectra of ir + , K + , and .ft' - . 
The data points come from the BRAHMS collaboration [32l . 
[3^ ] (circles - ir + , squares - K + , triangles - K~), while the 
histogram lines show the result of the model simulation with 
THERMINATOR. For n + the solid (dashed) line corresponds to 
the full feeding (no feeding) from the weak hyperon decays. 
The model parameters are from Eq. (|17|l . The experimental 
pion yields are corrected for weak decays as described in 



net protons. Since the p and p data carry no feed-down 
corrections for weak decays |32[ , one should compare the 
solid lines to the data. The shape of the p and p spectra is 
properly reproduced, but the model overshoots the data 
by about 50%. This feature occurs at all rapidities, also 
at mid-rapidity. The mismatch could be improved by 
decreasing T by a few percent and redoing the whole 
analysis, but we do not take the effort here, holding to the 
value (p~2|) obtained from global fits to all RHIC data for 
the particle yields at midrapidity. We provide, however, 
the results of the model calculation with no feeding from 
the hyperon decays, since it provides some measure of 
the systematic uncertainties in determining the proton 
and antiproton yields. 

Quite remarkably, the qualitative growing of the net- 
proton spectrum with y is obtained. This has a simple 
explanation on the ground of statistical models, since ap- 
proximately p — p ~ sinh(/xs(y)/T). Thus at RHIC the 
proton and antiproton spectra may be qualitatively ex- 
plained solely on the ground of the statistical approach. 

Figure [7] shows the pr-spectra at subsequent rapid- 
ity bins of the BRAHMS experiment. From top to bot- 
tom we have for pions y G [-0.1,0.0], [0.0, 0.1], [0.4, 0.6], 
[0.6,0.8], [0.8,1.0], [1.0,1.2], [1.2,1.4], [2.1,2.3], [2.4,2.6], 
[3.0,3.1], [3.1,3.2], [3.2,3.3], [3.3,3.4], [3.4,3.66], and 
for kaons y E [-0.1,0.0], [0.0, 0.1], [0.4, 0.6], [0.6,0.8], 
[0.8,1.0], [1.0,1.2], [2.0,2.2], [2.3,2.5], [2.9,3.0], [3.0,3.1], 
[3.1,3.2], and [3.2,3.4]. Each lower curve is subsequently 
divided by the factor of 2 in order to avoid overlap- 
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FIG. 6: (Color online) Top: the rapidity spectra of p and p. 
Bottom: spectrum of net protons, p — p. The data points 
come from the BRAHMS collaboration [32I . [3^ ]. while the 
solid (dashed) histogram lines show the result of the model 
simulation with THERMINATOR with full feeding (no feeding) 
from the weak hyperon decays. Data points should be com- 
pared to the model with full feeding (solid lines). The model 
parameters are from Eq. (|17[l . 



ping. The solid lines show the model calculation with 
optimum parameters (fT7|) . We note that the basic fea- 
tures of the experiment are reproduced, with the slope 
increasing with y. This can be explained with a lower 
transverse flow at larger y, as enforced by the parame- 
terization ([5]) . The quality of the agreement is similar in 
all rapidity bins. In Fig. [8] we give the similar study for 
the protons and antiprotons. The data points come from 
the BRAHMS collaboration and contain no weak- 
decay corrections, hence the solid lines should be com- 
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FIG. 7: (Color online) The p-r-spectra of ty + (top left), ir~ (top right), K + (bottom left), and K~ (bottom right) in the 
subsequent BRAHMS rapidity bins, see the text for details. The data points come from Ref. [3^], while the histogram lines 
show the result of the model simulation with THERMINATOR. 



pared to the data. Nevertheless, as in Fig. [SJ we also 
present the calculation with feed-down from the weak 
decays switched off, as it provides a measure of system- 
atic uncertainties. It should also be kept in mind that 
these uncertainties are quite large for the py-spectra of p 
and p, as can be inferred from the comparison of results 
of various experimental collaborations at RHIC (cf. for 
instance Fig. 12 of Ref. [H). 

At this point we have accomplished the goal of fixing 
the "fireball topography": we have the geometry/flow 
as well as thermal parameters dependent on the variable 
aii. Next, we may proceed as in the case of the boost- 
invariant model used at mid-rapidity, and compute many 
observables in addition to those already used up to fix 



the model parameters. These observables include one- 
body observables, such as spectra of various particles, 
including hyperons, mesonic resonances, etc., as well as 
two-body observables related to correlations: HBT radii, 
balance functions in rapidity, or event-by-event fluctu- 
ations. Here we only present a sample prediction for 
rapidity spectra of hyperons, shown in Fig. [9l An in- 
teresting feature is the very small splitting of f2 and 0, 
which results from the fact that /is — 3/is — 0, cf. Fig. [5] 



9 




P T [GeV/c] 



10 

E-, 

o, 1 

f lo 1 



10" 



lo 3 r 



10' 

















- " s "-o-,_o. 

3 T ~r 

* □ □ 


o 1 — 1 

'--i o 

□ 1 — 1 T 
x T T L Eh 


p 

""! o 










AAA 


—•-A 

A"'-A-^1 


A ■ p 


-P ( 


72; 






: P 


(/10) 


1 
























y~2 















E-H 

B, 

E-H 



10' E 



10- 



io- 4 , 



10 

S-H 

B, 















- 


' - - , 




P 






- 




gl - 










p (/20) 


i VTr 






- 






















, , , i 



0.5 



7.5 2 2.5 3 
P T [GeV/c] 



10 J i 







0.5 



1 



1.5 



Iff 















o 1 

--- | 
L„0 — 


^-1 1 p 

o\ -* 

■ : 








□ 1 — 1 
□ ' — 1 

--- □ 


□ ' — 1 i--7°~\ Q. 

~j, bsb j l., ]o 


6 








^Al 1_ ■ L mp _ 

— 1 


|_T 

s p. 










□ 

-ll- 

[A- 




\y»3 




1 




, , , , 1 







0.5 



1 



1.5 



2.5 3 



P T [GeV/c] 



2 2.5 3 
P T [GeV/c] 

FIG. 8: (Color online) The proton, antiproton, and net-proton transverse momentum spectra at y ~ 0, 1, 2, 3. The antiproton 
spectra and the net-proton spectra at y ~ 3 have been divided by 10. The data points come from the BRAHMS collaboration 
[33 ] and contain no weak-decay corrections. Solid (dashed) histogram lines show the result of the model simulation with full 
feeding (no feeding) from the weak hyperon decays. Solid lines should be compared to the data. 



IV. CONCLUSION 

The paper contains results of the single- freeze-out ther- 
mal model for rapidity-dependent spectra in relativistic 
heavy-ion collisions. We have used THERMINATOR to run 
the simulations and the BRAHMS data for ^Jsmn = 
200 GeV Au + Au collisions to fix the model parame- 
ters. Such a simulation is necessary when the system 
is not boost-invariant. It allows for an exact incorpora- 
tion of the space-time dependence of thermal parameters, 
precise inclusion of resonance decays, as well as incorpo- 
ration of experimental cuts. The extension of the original 
boost-invariant single-freeze-out model includes a modi- 



fication of the shape of the fireball, which here becomes 
narrower as the magnitude of the spatial rapidity an in- 
creases, as well as admits the dependence of the thermal 
parameters on a\\. As a result of a fit to the BRAHMS 
data we have obtained the dependence of the freeze-out 
chemical potentials on a\\. The freeze-out temperature 
is taken constant in the considered range of rapidities. 
With this extension we are able to properly describe the 
double d 2 N / (2 r KpT<lpTdy) spectra from the experiment. 
We also make predictions for other particles, in particular 
for hyperons. 

A code incorporating the elastic collisions neglected in 
the single-freeze-out approach could be used as an "af- 
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terburner" starting from our freeze-out condition. That 
way a more accurate collision picture could be achieved. 
As we have already mentioned, a recent study of Ref. [47| 
revealed that for the mid-rapidity p-r-spectra the elastic 
rescattering is not very important. 

Certainly, the scheme of this paper can be used for 
other collisions where departures from the boost invari- 
ance are significant, in particular for the rich SPS data. 
As we have said, the modeling involves the choice of the 
parameterization for the shape of the fireball and the ve- 
locity field of flow, where in fact we have quite a lot of 
freedom, as well as the dependence of the thermal param- 
eters at freeze-out on the space-time position. Accurate 
data for numerous observables as functions of the rapid- 
ity, not only abundances and spectra but also the correla- 
tion data (HBT radii, balance functions), would greatly 
help to constrain the freedom and acquire insight into 
the space-time evolution picture of boost-non-invariant 
systems formed in relativistic heavy-ion collisions. Most 
importantly, the knowledge of the dependence of i? S idc on 
y would put constraints on the shape of the fireball. 

Here are the main results of the paper: 

1. Naive extraction of the baryon and strange chemi- 
cal potentials from ratios ofp/p and K + / K~ works 
surprisingly well, as shown in the comparison to the 
full calculation in Fig. O 

2. The baryon and strange chemical potentials grow 
with a j_ , reaching at y ~ 3 values close to those of 
the highest SPS energies of ^Jsnn — 17 GeV. This 
agrees with the recent conclusions of Roehrich [3l| . 

3. At mid-rapidity the values of the chemical poten- 
tials are even lower than derived from the previous 
thermal fits to the data for \y\ < 1, with our val- 
ues taking fi B (0) = 19 MeV and [i s (0) = 5 MeV. 



The reason for this effect is that the particle with 
\y\ < 1 originate from a region |aii | < 2, and on the 
average the effective values of chemical potentials 
are larger compared to the values at the very origin 
(cf. Fig.©. 

4. The local strangeness density of the fireball is com- 
patible with zero at all values of an. Although 
this feature is natural in particle production mecha- 
nisms, here it has been obtained independently just 
from fitting the chemical potentials to data. 

5. The ratio of the baryon to strange chemical poten- 
tials varies very weakly with rapidity, ranging from 
~ 4 at midrapidity to ~ 3.5 at larger rapidities. 

6. The d 2 N/ (2Trp±dp±dy) spectra of pions and kaons 
are well reproduced, supporting our hypothesis for 
the shape of the fireball in the longitudinal direc- 
tion. 

7. The rapidity shape of the spectra of protons and an- 
tiprotons measured by BRAHMS [32| is described 
properly, while the model predict too large nor- 
malization, overproducing these particles by about 
50%. This suggests a lower value of T by a few 
percent, or presence of non-equilibrium factors. We 
also note that the feature of an increasing yield of 
the net protons with rapidity is obtained naturally, 
explaining the shape of the rapidity dependence on 
purely statistical grounds. 
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